Dynamical enhancement of spatial entanglement in massive particles 
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We discuss dynamical enhancement of entanglement in a driven Bose-Hubbard model and find 
an enhancement of two orders of magnitude which is robust against fluctuations in experimental 
parameters. 
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Quantum coherence is often thought to be found only 
in very small systems or under artificial lab-conditions, 
since otherwise unavoidable environment coupling results 
in rapid loss of coherence. Whereas this holds quite gen- 
erally, it holds in particular for many-body coherence, 
synonymous for entanglement. As recent experimental 
[H and theoretical [H evidence suggests, however, excep- 
tions to this general rule exist. In particular coherent 
driving can compensate for the environment-induced loss 
of coherence and, thereby stabilize entanglement under 
conditions under which a static system would be com- 
pletely separable 0, HJ ■ 

Such dynamically induced entanglement does not only 
hold the potential to influence macroscopically observ- 
able properties [5j , but certainly also opens up new paths 
towards scalable quantum information processing which 
otherwise is limited through the unfavourably scaling de- 
phasing times with the system size @, 0|- However, our 
current understanding of dynamical enhancement of en- 
tanglement is still in its infancy. 

In this Letter, we investigate the dynamical enhance- 
ment of entanglement between ultracold bosonic atoms 
stored in an optical lattice that gives rise to a spatially pe- 
riodic potential created by two counterpropagating laser 
beams of wavelength A and amplitute Vq in one direction. 
A tight perpendicular confinement of strength V±_ in the 
other two directions restricts the motion of the atoms to 
one dimension. In the deep lattice limit Vb 3> En, where 
Er = h 2 k 2 /2m (with k = 2ir/X) is the recoil energy, and 
at sufficiently low temperatures, this system can be well 
described l8fl in terms of the Bose-Hubbard Hamiltonian 



L-l 



H 



k+i + a] +1 ai) + 



U 



1=1 



L 

E 

i=i 



ni("i-l)(l) 



where the creation operator aj creates and the annihi- 
lation operator a; annihilates a boson at lattice site 
The tunneling parameter J and the on-site interaction U 
depend on the lattice parameters approximately [9( via 
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Whereas the lattice depth Vb is typically time- 
independent, we will compare here the dynamics of such 
an autonomous system with its driven version (Io| , where 
Vb is modulated temporally 



Vb(t) = v(l + dVsin(wi)) . 



(4) 



As we will show, the temporal modulation of the tunnel- 
ing parameter J and the on-site interaction U that results 
from this lattice depth modulation drives the atoms into 
a spatially strongly correlated, i.e. entangled state. For 
the verification of the entanglement properties, we envi- 
sion a rapid separation of the few-body system into two 
parts, what can be realized by ramping up a potential 
barrier as depicted in Fig. [TJ 
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FIG. 1: (color online) a) The system is prepared in an en- 
tangled state, b) The system is split into two halves by raising 
the intermediate barrier. J is the kinetic term and U the on- 
site interaction in the Bose-Hubbard model |T]l. 

This spatial separation effectively switches off the in- 
teraction between the two subsystems what freezes the 
entanglement dynamics. However, it also entails that 
each subsystem will typically not have a well-defined 
particle number. Whereas correlations in the particle 
number formally may imply entanglement, its experi- 
mental verification will be technically impossible because 
this would require the measurement of coherent super- 
positions of states with different numbers of massive par- 
ticles [TlJ. Therefore, we will postselect 12 1 states with 
well-defined local particle number. We will focus in par- 
ticular on those cases in which particles are split evenly 
between the two subsystems, since this is the case that 
occurs with highest probability, and this is also the case 
in which the highest entanglement can be achieved. Do- 
ing so, we obtain a clean notion of entanglement between 
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the two separated halves of the optical lattice where each 
half is filled with a fixed number of particles and can be 
addressed individually. 

In the following we quantify the entanglement of the 
ostselected states with the entropy of entanglement 



1 31 ] in the case of pure states and the negativity [14j 
in the case of mixed states. The entropy of entangle- 
ment is given by the von Neumann entropy of the re- 
duced density matrix g r that is obtained through the 
partial trace over one subsystem of the entire many- 
body state, i.e. E(^) = — Trg r log 2 g r . The negativity 
N(p) = (||p PT || — l)/2 of a mixed state g is defined in 
terms of the trace norm of the partially transposited den- 
sity matrix p 
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FIG. 2: (color online) a) Entropy of entanglement E and 
probability for successful postselection of evenly distributed 
number of bosons between the right and the left part of a 
L = 6 well Bose-Hubbard system. The overall particle num- 
ber is TV = 6 (solid blue), TV = 3 (dashed red) and TV = 2 
(dash-dotted green), b) Negativity for the thermal state at 
temperature T — nK (solid blue) , T = 40 nK (dashed red) 
and T = 80 nK (dash-dotted green) for TV = 6 particles. 

Fig-Hk) depicts the ground state entanglement proper- 
ties for the exemplary case of L — 6 wells filled with 
TV = 6, 3 and 2 bosons, respectively. In all these cases 
the system features extremely low entanglement, which is 
also observed for different system sizes and particle num- 
bers. The qualitative dependence of E on the parameter 
U I J is as expected: If the tunneling dominates the sys- 
tem dynamics, i.e. if U/J ~ 0, the bosons populate the 
same single-particle states such that after postselection of 
local particle number the system is separable. For finite 
interaction U the bosons repel each other and establish 
correlations; therefore, E typically increases with U/J. 
There are, however, exceptions, like the case of unit fill- 
ing (depicted in blue), where a separable, perfect Mott 
insulator 0, ||| develops for U ' / 'J — > oo. This is also re- 
flected in the fact that in this limit the bosons will always 
be separated in a balanced fashion between the left and 
right half of the system, whereas typically the probability 



for this is smaller than unity. 

Fig. [2] shows data for rather small fillings, but also the 
entanglement properties of the ground state for larger 
filling N/L > 1 can be inferred from this data since sys- 
tems with TV + mL bosons (with integer m) behave qual- 
itatively similar as the system with TV bosons. That is, 
the maximal amount of entanglement in the ground state 
does typically not exceed the value of E « 0.05. 

Assuming perfect ground state cooling is certainly 
a theoretical idealization, but also thermal excita- 
tion cannot enhance the entanglement as shown in 
Fig. [^b), where the negativity of the thermal state 

pth = ex P(j^pf)/Z i s shown. The probability to find the 
bosons split evenly into left and right half decreases with 
increasing temperature, and the entanglement is even 
lower than for T = 0. 

As we will see in the following, this is strongly con- 
trasted by the behavior of the driven system. To 
be specific, we consider the experimental parameters: 
V = 10 E R and V±_ = 30 E R as lattice depths, dV = 0.2 as 
lattice depth modulation, A = 842 nm as the wave length 
of the laser, a s = 5.45 nm as the scattering length and 
m = 86.909 u as the mass of rubidium-87 [15]. As driv- 
ing frequency, we chose w = U/h = 12862 Hz, what corre- 
sponds to resonant driving in the Mott-insulating regime. 
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FIG. 3: (color online) Entropy of entanglement E and proba- 
bility of successful postselection for the driven Bose-Hubbard 
Hamiltonian © with L = 8 = TV (solid blue), L = 6 = TV 
(dashed red), L — 4 = TV (dash-dotted green), L = 6 and 
TV = 5 (thin dash-double-dotted orange), and L = 6 and 
TV = 7 (thin dotted turquoise) . At t = 100 ms the driving is 
stopped and at t = 150 ms the lattice depth is increased in 
order to freeze the entanglement dynamics. 

Fig. [3] shows the dynamical enhancement of entangle- 
ment caused by the driving. We start with the ground 
state of the static system. After 100 ms the driving is 
switched off and after t = 150 ms the lattice depth is in- 
creased to Vo = 30 Er in order to freeze the entanglement 
dynamics completely. Apparently, entanglement grows 
rather quickly once the driving is switched on. After 
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about 10 to 20 ms the increase slows down a bit, until 
entanglement saturates after t w 50 ms. From that time 
on, entanglement fluctuates around an average value due 
to the finite interactions in the system. Fig. [3] shows 
how these remaining fluctuations smooth out with grow- 
ing system size. Whereas the saturation time turns out 
to depend crucially on the lattice depth modulation dV, 
the final values of entanglement seem not to. Thus, a 
larger modulation of the lattice depth can accelerate the 
entanglement generation; however, even very weak driv- 
ing can yield the same enhancement of entanglement as 
strong driving. 

Once entanglement is saturated, it can be preserved 
with the help of a potential barrier that separates the 
two entangled atomic ensembles as depicted in Fig. Q] 
and effectively switches off their mutual interaction. The 
small fluctuations around the average entanglement that 
persist for t > 100 ms can be stopped with an increase of 
the lattice depth, here at t = 150 ms. 

As it can be seen in Fig. [3^) the attainable entangle- 
ment grows with increasing system size. This is as ex- 
pected since larger systems can carry more entanglement. 
A concern, however, is that with increasing number of 
atoms, also the number of possible distributions of atoms 
between the two subsystems is growing, so that the prob- 
ability to find an even distribution might decrease. This 
decrease is apparent in Fig. 13b) , where the probability of 
an even distribution drops from p — 0.33 for 4 atoms to 
p = 0.27 for 5, 7 and 8 atoms. It also becomes apparent, 
however, that this decrease occurs in small systems, and 
that these probabilities become largely independent of 
the particle number for large N. Thus, one should expect 
to find evenly distributed particle numbers with substan- 
tial probability also in an experiment with significantly 
more bosons than a numerical simulation can handle. 

In an experiment, certainly also the timing will be cru- 
cial. As Fig. [3] shows, the system evolves rapidly, and 
fluctuations in the durations of driving or ramping up the 
barrier that are comparable to system time scales will re- 
sult in the generation of a mixed state which typically has 
reduced entanglement. The relevant time scale can be 
obtained from the fidelity f(t, At) = \(ip(t)\ip(t + Af))| 2 , 
where \if>(t)) is the postselected system state after driv- 
ing of duration t. The fidelity / is depicted in Fig. 0] 
for the exemplary case of N — 6 particles in a L = 6 
site lattice. The width of the central peak (full width 
at half maximum) that determines the minimal required 
experimental precision reads in this case At m = 0.1ms. 
In a similar fashion, we can also estimate the required 
precision for all other experimental parameters, such as 
the potential V (AV m = 0.08 Er), the perpendicular con- 
finement V± (AVj_ )TO =0.12Er), the amplitude of the 
driving dV {AdV m = 0.016) and the driving frequency u) 
(Aw m = 14.5 Hz). 

To estimate the impact of fluctuations of these para- 
meters on the attainable entanglement, we have to con- 
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FIG. 4: a) Fidelity f(t, At) at t = 100 ms for state preparation 
with imperfect timing, b) Negativity N of the mixed state ([5]) 
resulting from the measurement of the driven Bose-Hubbard 
Hamiltonian with L — 6 = TV at to = 100 ms. 

sider the mixed state that is obtained with many repe- 
titions of the experiment with the fluctuating parameter 
taking different values at each repetition. To be specific, 
we focus on the inaccuracy in the duration of driving, 
and we assume that these durations are distributed ac- 
cording to a Gaussian centered around to = 100 ms with 
a standard deviation r. This gives rise to the mixed state 

p(r) = — L= f°° dt e-^m))(m\ • (5) 

The negativity of this mixed state is depicted in Fig. HJd) 
as function of the inaccuracy r of the duration of driving. 
For T=0ms, the situation reduces to the case of pure 
state entanglement as discussed above; but for finite r en- 
tanglement is reduced significantly what is a very generic 
feature of mixed states. Besides this expected behavior, 
there are two features that should be stressed: 

First, at r = ms the first derivative of N(t) van- 
ishes, so that entanglement turns out to be insensi- 
tive to small timing errors. The second-order Tay- 
lor expansion reads N(t) « 3.064 — 900(r/ms) 2 . 
That is, timing errors below 0.01 ms imply a change 
in negativity of less than 3 %. 

Second, even in the presence of significantly larger 
timing errors there is still rather strong entangle- 
ment with N(t) ~ 1, i.e. a value attained for a 
maximally entangled Bell pair! 

In particular, this astonishing robustness against experi- 
mental fluctuations underpins that potential that driving 
offers as means to create entanglement as compared to 
engineered interactions. 

As recent investigations on driven spin-systems suggest 
, the feature of dynamical enhancement of entangle- 
ment is not particular for the Bose-Hubbard system, but 
a rather generic feature, that is largely independent of 
detailed system properties. An advantage of the present 
bosonic system as compared to many spin-systems is that 
particle numbers can easily be varied in an experiment, 
what provides the means to study the generation of en- 
tanglement in the entire regime from rather small systems 
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through the mesoscopic domain, up to the semi-classical 
regime. In particular, observing the rise and decay of en- 
tanglement with increasing particle number will provide 
us with valuable insight in the emergence of classical be- 
havior in large quantum systems. 
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